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Abstract 



The NRQCD factorization formalism of Bodwin, Braaten, and Lepage pre- 



scribes how to write quarkonium production rates as a sum of products of 



short-distance coefficients times non-perturbative long-distance NRQCD ma- 



^ . trix elements. We present, in the true spirit of the factorization formalism, a 

detailed calculation of the inclusive cross section for hadronic ip production. 
We find that in addition to the well known color-singlet production mecha- 
nisms, there are equally important mechanisms in which the cc pair that forms 
the ip is initially produced in a color-octet state, in either a ^5*1, ^5*0, "^Pq or 
angular-momentum configuration. In our presentation, we emphasize the 
"matching" procedure, which allows us to determine the short-distance coef- 
ficients appearing in the factorization formula. We also point out how one 
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may systematically include relativistic corrections in these calculations. 



Typeset using REVT^iX 



I. INTRODUCTION 



Bodwin, Braaten, and Lepage (BBL) have devised a rigorous factorization formalism 
that places the calculation of inclusive quarkonium annihilation and production on a solid 
theoretical foundation [jl|. Their approach is based on nonrelativistic quantum chromody- 
namics (NRQCD) 0. This is an effective field theory involving an expansion in derivatives, 
constructed so as to be equivalent to QCD for non-relativistic heavy quark-anti-quark scat- 
tering, to any desired order in the relative three-momentum of the heavy quarks. In Ref . [|T| , 
it is shown that the production rate for a quarkonium bound state H can be written using 
a factorization formula, which is a sum of products, having the form 

^(^) = E^(o|c?f|o). 

The "short-distance" coefficients F„ are obtainable in perturbation theory as a series in 
as{mQ), and can be determined through a calculation of the production rate of an on-shell 
QQ pair. As to the NRQCD matrix elements (0|(9^|0), they encode analytically uncalculable 
"long-distance" effects such as the hadronization of a QQ pair into bound quarkonium. The 
index n labels color and angular-momentum quantum numbers, as well as the order in the 
NRQCD momentum expansion. In contrast to previous approaches, the BBL formalism 
establishes a framework within which it is possible to properly handle soft gluon effects and 
systematically incorporate relativistic corrections. 

The size of any NRQCD matrix element can be estimated as some power of the small 
parameter f , the typical velocity of the heavy quarks in the bound state. (For charmonium, 
f ^ ~ 0.3.) Combining f-scaling estimates with knowledge of the size of the Fn (for which the 
order in as(mQ) is easily determined), one can deduce the relative importance of the various 
terms in the factorization formula. Keeping only those terms warranted by experimental 
precision, one can cast each observable as a sum involving only a small number of NRQCD 
matrix elements. In this way, the important NRQCD matrix elements can be empirically 
determined by fitting a set of factorization formula to a body of experimental data. 
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In this paper we carefully present an example of a calculation in the factorization for- 
malism. The quantity that we choose to calculate is the inclusive hadronic ip production 
cross section, since it is needed for a complete theoretical description of ip production at 
fixed target experiments. Using f-scaling rules alluded to above, we find that in addition to 
previously calculated color-singlet contributions, there exist equally important contributions 
involving the production of cc pairs in a color-octet state. The color-octet processes which 
turn out to be important in hadronic ip production are those in which the heavy quark pair 
is produced in either a ^Si, ^Sq, ^Pq or angular-momentum state. 

Taking the octet-mechanism contributions to hadronic production as an illustrative 
example, we present a systematic procedure for determining the short- distance coefficients 
Fn appearing in the factorization formula. This procedure involves the calculation of rates 
for the production of QQ pairs with specific color and angular-momentum quantum numbers. 
These production rates are expressed as an integral over q, the relative three- momentum of 
the Q and Q in the QQ rest frame, with the integrand being a Taylor expansion in q^/m^. 
The production rates are calculated both in full QCD and in NRQCD, and then the results 
are "matched," allowing a determination of the F„. 

The Feynman diagrams for the leading color-singlet hadronic ip production subprocesses 
are shown in Figure 1. However, we do not calculate these leading color-singlet contribu- 
tions since they have already been treated using the color-singlet- wavefunction model 0. 
For calculations at leading order in v"^, the results of this latter approach can be readily 
transformed into the factorization formalism language. 

II. NRQCD FACTORIZATION APPROACH 

A. Factorization Formula 

As was stated in the introduction, the NRQCD factorization formalism prescribes a 
means of expressing the inclusive production rate of a heavy quarkonium meson H using a 
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factorization formula, which is a sum of products of the form 

a{A + B^H + X)=Y: ^ (0|Of |0) . (2.1) 

n 

The short-distance coefficients F„ can be calculated using Feynman diagram methods. Each 
Fn can be expressed as a perturbation series in asirno). The appropriate scale for the 
process is irtq since this is the scale that is associated with the production of a QQ pair with 
small relative momentum. 

The NRQCD matrix elements (0|Cf |0) are of the form 

(0|O^|0) = EE mn\H^j+X) {H^^+X\1C'M , (2.2) 

X rnj 

where /C„ and /C^ are bilinear heavy-quark field operators; /C„ may take forms such as xV; 
xV*'?/', x^(^' D "^A, e-tc. and /C^ takes forms such as V'^X; cic. The symbol dn appearing in 
the factorization formula is the combined mass dimensions of the operators /C„ and K,'^. The 
heavy-quark field operators ip and ^ire color-triplet columns, and are defined to annihilate 
a heavy quark and anti-quark, respectively. Thus, the /C„ annihilate QQ pairs while the /C^ 
create them. The insertion of a color matrix T*^ in the bilinear operators /C„ and /C^, e.g. 
X^T°"ip, means that the bilinears project a color-octet rather than a color-singlet QQ state. 
The index n labels the various properties of the /C„ and /C^. These properties include: 1) 
the color-state of the QQ pair projected by the operator (singlet or octet, designated by 
1 or 8); 2) the angular- momentum state of the QQ pair projected by the operator (given 
using the spectroscopic notation "^^^^Lj); and 3) the order of the operator in the momentum 
expansion in the NRQCD effective lagrangian. The order in the momentum expansion can 
be increased simply by inserting into /C„ or /C^ the scalar (— | T))'^/ml, where D is an 
SU{3) covariant derivative. In Eq. ( p.2| ), the sum over X is a symbolic reminder that we are 
calculating inclusive quarkonium production. In the discussion that follows, the sums over 
X and mj will always be assumed. 

The factorization formula Eq. ( |2.1| ) contains an arbitarary factorization scale A. In the 
matrix elements the scale A can be identified with the ultraviolet cutoff of the NRQCD 



effective theory. Since physical results are independent of A, dependence on A in the short- 
distance coefficients cancels against that in the NRQCD matrix elements. 

In the calculation of the factorization formula for a given quarkonium production rate, 
it is first necessary to determine which are the most important terms in the series. This 
operation entails a consideration of two issues: the order in as{mQ) of the F„ and the order 
in of the (0|(9„|0), i.e. their "f -scaling." In the next subsection, we will turn to the issue 
of f-scaling of the NRQCD matrix elements. 

B. w-scaling of NRQCD matrix elements 

The most practical approach to the issue of f-scaling of NRQCD matrix elements is 
to define a f-scaling "baseline," with all other matrix elements being suppressed by some 
relative power of v. This baseline corresponds to those matrix elements in which the bilinear 
operators /C„ and /C^ project out the predominant QQ component of the quarkonium, and 
in which this predominant component is in an S-wave. For charmonium, examples of such 
baseline matrix elements are 

(0|O?(%)|0) = (0|xVl^ + X) {v + X\ij^x\0) 

{0\OtCSi)\0) = {0\x^cT^\^Pm,+X) ■ (V^™,, +X|^t^^|o) . (2.3) 

If one inserts the scalar operator (— | D)^/ ml in either of the bilinear operators, the resulting 
matrix elements are suppressed by compared to those in Eqs. ( p.3|) . 

Next let us consider those matrix elements in which the bilinear operators project out the 
predominant QQ component of the quarkonium, and in which this predominant component 
is in a P-wave. Examples are 

{0\OlCP,)\0) = (0|xt(-^D)vlV.+^)-(V.+^l^^(-^D)x|0) (2.4) 
(0|OfePo)|0) = ^(0|x^(-^D-^)^|xo + X)(xo + X|^t(^_^5.^^^|0) (2.5) 

{0\OrCPim = ^(0|x^(-^Dxa)^|xi™,+X)-(xi^,+X|^t(^_^5x^)x|0) (2.6) 

6 



{0\OrCP2m = {0\x^ {-'- a^-l) ij\x2n.j+X) {X2mj+X\ij^ (-^ a^-l) x|0) • (2.7) 

To a rough approximation, the above matrix elements, when divided by ml, are suppressed 
with respect to the basehne by v'^, and of course, the insertion of (— | D)^/m^ would result 
in further suppression by v'^. 

The preceding discussion has merely involved counting powers of momentum. There is, 
however, an important additional issue in the v-scaling of NRQCD matrix elements. This 
second issue hinges on the Fock state expansion of the quarkonium meson. The Fock state 
expansion of, for example, the ip particle, in Coulomb gauge, can be thought of schematically 
as 

I V' ) = ^c-|cc(l, ^Si) ) + A,sg\cc{8, ^Pj) g) + A,-,gg\cc{S, ''S^) gg) + B,^gg\cc{l, ^S^) gg) 

+Cccgg|cc(8, ^Dj) gg) + D,^gg\cc{l, ^Dj) gg) + B,,g\cc{8, ^So) g) + ■ ■ ■ , (2.8) 

where g represents a dynamical gluon, i.e. one whose effects cannot be incorporated into 
an instantaneous potential and whose typical momentum is rricv'^. The angular-momentum 
quantum numbers of the cc pairs in the various Fock states are indicated in spectroscopic 
notation, and their color configurations are labeled by 1 for singlet or 8 for octet. We will 
now discuss the v-scaling of the various coefficients A, B, etc. 

Since the state \cc{l,^S\)) is the predominant Fock state in 1-0), we expect that the 
coefficient Acc is just a little less than unity, i.e. Acc ~ v^. As to the Fock state |cc(8, ^Pj)g), 
this configuration arises when the predominant state radiates a soft dynamical gluon; such 
a process is mediated principally by the electric dipole operator, for which the selection rule 
is L' — L ±1, S' — S , and which involves a single power of heavy quark three- momentum; 
thus, the coefficient Accg is of order v^. The electric dipole emission of yet another gluon 
involves a change from the P-wave state Accg\cc{^,^Pj)g) to the S- and D-wave states 
\cc{%,^Si)gg), \cc{l,^Si) gg), \cc{S,^Dj) gg) and \cc{l,^Dj) gg); their coefficients — A^cgg, 
Bccgg, Cccgg and Dccgg — are of order v^. Lastly, we consider the coefficient of the state 
\cc{8,^SQ)g). Fluctuations into this spin-singlet state from the predominant spin-triplet 
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state are associated with the emission of a soft gluon via a spin-flipping magnetic dipole 
transition; such transitions involve the gluon three-momentum (~ rricV^) rather than the 
heavy quark three-momentum (~ rricV^), and therefore the coefficient Bccg is of order v^. 

We now bring together the power-counting rules and the Fock state ideas presented 
above. Let us consider as an example the matrix element 

(oiOg^f^Olo) = (o|xVr«V'|V'rn, + ^)(V'm, + ^IV'Vr'^x|o). (2.9) 

The bilinear operators project out the \cc{S^^Si)gg) Fock state, whose coefficient in the 
expansion is of order v"^. Thus, (OjOg (^5'i)|0) is suppressed by roughly v'^ compared to the 
5'-wave baseline. 

As another example, consider the v-scaling of the matrix element {Q\0^{^ Pj)\Qi) . The 
bilinear operators contained therein project out the state |cc(8, ^Pj)gi), whose coefficient 
in the expansion would be of order v^. This Fock state suppression, combined with the 
derivatives inherent in the P-wave operators, gives a v-scaling of (with respect to the 
baseline) for (0|e>^(3Pj)|0). 

The v-scaling of those NRQCD matrix elements pertinent to the calculation of hadronic 
■0 production are given in Table I. 

C. Order in ag of short-distance coefficients 

The leading-order Feynman diagrams for the production of a cc pair in color-octet states 
arc presented in Figures 2 and 3. The diagram in Figure 2 [qq g ^ cc) produces 
a cc(8, ^Si) state. The diagrams in Figure 3a (the gluon fusion process gg — > cc) produce 
cc(8, ^Sq), cc(8, ^Si) and cc(8, ^Pj) states. The diagram in Figure 3b {gg g ^ cc) produces 

a cc(8, ^Si) state. 

It is interesting to note that, at leading order in the q^/m^ expansion, the amplitude 
for the production of cc(8, ^5"!) from Figures 3a cancels against the amplitude for the same 
quantum number production in Figure 3b, so that the total production of cc(8, ^Si) from 
gluon-gluon collisions vanishes. 

8 



To determine the rough size of a term in the factorization formula, one attributes to it a 
power of as^rric) ~ 0.3 for each vertex in the Feynman diagram, and the appropriate power 
of 0.3 according to its v-scahng. We summarize the as{mc) powers associated with the 
leading order subprocesses in Table II. 

Once the most important terms in the factorization formula have been identified, one 
procedes to compute the short-distance coefficients F„. This is done by matching pertur- 
bative full QCD and perturbative NRQCD calculations for the production rate of cc pairs 
with specific color and angular-momentum quantum numbers. These rates are expressed as 
an integral over relative momentum q, with the integrand given as a Taylor expansion in 

In the next sections, we calculate the production rates for the various important octet 
subprocesses contributing to hadronic ip production, and outline in detail the matching 
procedure. 



III. PRODUCTION RATE FOR qq ^ i; SUBPROCESS 

We now calculate the leading-order contribution to the cross section for inclusive ip 
production from the color-octet subprocess qq — > ■0. This is done in three steps. First we 
carry out a perturbative QCD calculation of the production rate for the process q{k)q{k) 
c{p)c{p) with on-shell heavy quarks. The second step consists of calculating the production 
rate for the same process in NRQCD. The final step will be to determine the short-distance 
coefficient F„ by matching the QCD and NRQCD calculations. 



A. Production rate in full QCD 

In noncovariant conventions, the transition amplitude for the process q{k)q(k) — > 
c(p)c(p), illustrated in Figure 2, is given by 

r = y mrTMk) Usiph.T^^vM , (3.1) 

4.^kokoPoPo [P+P) 
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where I, I, m and m are color indices and s and s are heavy-quark spin indices. The matrices 
are normahzed throughout so that Tr[T°'T''] = 6°'^/2. It is convenient to re-express the 
c and c momentum in terms of q (the relative three-momentum of the quark and antiquark 
in the cc rest frame), and P (the total four- momentum of the quark and anti-quark in the 
lab frame): 

= ^P^ + K^\P)q\ 

= ^P" - K^'\P)q\ (3.2) 
where the Lorentz boost matrix A^'^(P) is given by 

AO = J_p3 

pi pj pO pi pj 

with E"^ = rn^ + q^. The boost matrix has the following useful properties: 

pfj. pu 



A'^A'^g,, = -6,, . (3.4) 

The heavy quark Dirac bilinear u{p)'y'^v{p) appearing in Eq. ( p.l|) is now re-expressed in 
terms of q, the Pauli two spinor ^ and antispinor rj, and the Pauli matrices a* [^: 

Us{p)rv,{p) = a; (2E, CV^r/, - q'^k " ^v)j , (3-5) 

where Us{p)Us{p) = 2mc and ^l^t = ^st- The procedure of re-expressing the manifestly 
Lorentz-invariant Dirac bilinear current in terms of two-spinors can be called the "reduction 
of Dirac bilinears." Applying the reduction to Eq. (|3.1| ) we obtain 



2v/A:o^o(Po'-4(Ao-q)2) 



(3.6) 



(In the above context, the symbol represents a number, not a matrix; in our notation, 
the two-spinors ^ and 77 do not carry a color index and are not acted upon by the T"" . 
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Our strategy of placing the factor between the two-spinors is designed to make the 
matching procedure more transparent.) The bihnear combination C,l<yji]s corresponds to a 
^5*1 spin configuration of the cc pair. Working to leading order in q, we keep only the ^Si 
term: 

-r = / - , v{k)^^Tf,u{k) ^ A,,(P) CV^T^^r/, . (3.7) 

Anticipating that the probability of hadronization is negligible for |q| > ^qcd, we need only 
derive the algebraic form of T for small q, and therefore need only keep the first few terms 
in the Taylor expansion in q^/m^: 

T ^ ^l— vCkh'T-tMk) ^ A,, CV^T^^r^s- f 1 + a4 + ■ ■ ■) , (3-8) 

where the constant a — which we do not bother to calculate explicitly — is the coefficient 
of the second term in the Taylor series expansion in q^/m^. It must be kept in mind that 
the factor appearing in brackets is not, at this point in the calculation, a rapidly converging 
series. Indeed, the upper bound on q is associated with the total momentum available in the 
experiment, and in general this can be much greater than rric. However, for the purposes of 
matching in the factorization formalism, we treat T as an expansion in q^/m^, and consider 
only the first few terms, the probability of hadronization being non-negligible only for small 

q 

Multiplying ( |3.8| ) by its complex conjugate, summing over final colors and spins, and 
averaging over initial colors and spins, we obtain 

^ ■ vWT;^m^s (l + 2a^ + . . , (3.9) 

where a sum over spins s and s and color m and m is assumed. Here we have made use of 
the fact that 

ss ss 

The expression for the cross section in terms of the transition amplitude squared is 
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where F = 2 is the noncovariant flux. Changing variables from p and p to P and q, the 
above expression becomes 

= TT 



"^'"^ 5iEf-E,)\r\' fl + «'^ + ---V (3-12) 



F J (27r)3 J (27r)3 2Eg 

(27r)3 "^"^ "'^ ' ' ^ m 
Inserting Eq. ( ^.91 ) into Eq. ( p.l2| ) we obtain the rate for on-shell cc production: 

a{qq - cc(8, ^^0) = ^ / ^(^/ " ^^.) ^^T^^^. V^Tg^^s (l + a"^ + ■■■). 

(3.13) 



B. Production rate in NRQCD 

The next step is to calculate the production rate for the process qq cc(8, ^S'l) in 
perturbative NRQCD. The flrst two terms are 

a{qq-^cc{8,'S^))= f-^( {0\Of^\'S,m + {o\vf^\'s,m + -- 

(3.14) 

where ^8(35*1) and G8{^Si) are speciflc cases of the short- distance coefficients and where 
{0\Of'\'S,)\0)= Y: (0|xVT>|c(.,m,q)c(.-m,-q)) 

ssmm 

{c{s, m, q)c(s, m, -q) \^^a'T''x\0) (3-15) 
{0\Vf'\'S,)\0) = Iy: (OIxVT^ (-^ D)%|c(.,m,q)c(^,m,-q)) 

{c{s, m, q)c(s, m, -q) |V^V'T"x|0) 
+ ^ E (0|xV'T>|c(s,m,q)c(s,m,-q)) 

ssmm 

(c(s, m, q)c(s, m, -q)|^VT'^ (-^ d)' x|0) ■ (3.16) 
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Our task here is to derive explicit expressions for the cc(g) matrix elements in Eqs. (p.l5|) 



and ( p.l6| ). For practical reasons, it is worthwhile to quickly note one possible set of con- 
ventions for these calculations. Let us define the single-particle annihilation and creation 
operators to obey the anticommutation relation [a(s, z, q), a^(t, j, p)] = (27r)'^5(q — p)5'^*5*-'. 
The annihilation operator acts according to a{s,i,q)\c{t,j,p)) = (27r)^(5(q — p)(5^*5*-'|0), 
and the creation operator according to a"''(t, j, q)|0) = |c(t,j, q)). Single particle states are 
normalized so that (c(s, q)|c(t, j, p)) = (27r)'^5(q — p)6^^6^K Then, conceiving of the field 
operators ip and x color-triplet vectors, we write their Fourier decompositions as 



(2vr) 

XU^) =T.j ^K^-, ^, P)^|e-^^-^ , (3.17) 

where and 77! are the two-spinors. 

Using the conventions outlined in the preceding paragraph, we derive explicit expressions 
for the right-hand-side of Eq. (|3.14|) , obtaining 



aiqq ^ CC(8, ^1)) = J ( + + ' ' ' ) ^>^rnnrVs VWT^m^s ■ 

(3.18) 



C. Matching 

Finally we determine the short-distance coefficient by matching the perturbative QCD 
result Eq. (|3.13|) to the perturbative NRQCD result Eq. 



2«27r3 

i^8(=^5i) = ^5(^-4m^) (3.19) 

G8(^^i) = ^^^5(^-4mD. (3.20) 

The matching procedure for on-shell scattering fixes the short-distance coefficient in the 
non-relativistic effective theory. However, these same coefficients apply to the calculation 
of the production of bound quarkonium states. This is because the non-perturbative effects 
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that will bind the c and c into a charmonium meson take place over much longer distances 
than the separation l/irtc associated with the formation of the cc pair. Thus, although 
the short-distance coefficient has been determined using a perturbative calculation of the 
production of free quarks and antiquarks, the same short-distance coefficient applies to the 
formation of a cc pair in a charmonium state. The nonperturbative effects involved in the 
formation of the boundstate are described by the matrix elements. Therefore, the inclusive 
cross section for qq ^ ip via the color-octet mechanism is, to leading order in f ^, 



a{qq^^) = ^^^^^{m('Sm 
ml 



5{s-AmlmotCsm- (3.21) 



27m3 

This procedure, which involves the determination of the -F„ for the production of on-shell 
heavy quarks, correctly takes into account the effects of binding energy, as can be verified 
by explicit calculations in NRQED j^. 

IV. CROSS SECTION FOR gg ^ ^p SUBPROCESS 

Next we turn our attention to the process g{gi)g{g2) c{p)c{p) which, at leading order 
in proceeds through the Feynman diagrams in Figure 3. To determine the short- distance 
coefficients we follow a sequence of steps analogous to that outlined in the previous section. 
At this point we will abandon or habit of writing heavy quark spin and color indices explicitly. 

The transition amplitude calculated from the diagrams in Figure 3a is 

T{gg cc) = _ -. e^(^i)e^(^2) 

4Vfl'io5'2oPoPo 



u{p) 



yayfe 7^(:^- ^l+^ch" j^bj^a Yi- P+ ^l+mc)r . . . ,^ 



The product of color matrices in Eq. ( |4. 1| ) can be rewritten as 



j.aj.b ^ l^ab ^ l^^afcc ^ ifabc^T", (4.2) 

6 2 
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where d"''"^ is a symmetric tensor and f"''"^ is an antisymmetric tensor. Since we wish to 
determine the short- distance coefficient for color-octet production, we will discard the 5"''/6 
term of Eq. (|4.2|). Then, using the Dirac equation, Eq. (|4.1|) can be transformed into 



T{gg cc) 



{91 ■ 



El , 



(^i,{.9i)(^u{92) 



_ 91 ■ Ag .^,u 
El ^ . 



where 



fiu _ ■ tJ.ua 



ie^'''''{gi-g2)alpl, + 2 [(Ag)^'^ + (Ag) V] , 



and 



We reduce the Dirac bilinears using Eq. ( p.5| ) and 



rric 



>v{p) , 
(4.3) 



(4.4) 



(4.5) 



(4.6) 



Keeping only terms of leading order in q, the resulting amplitude may be expressed in the 
form 



T{gg cc) ^ Tcs{gg cc) + Tcas{99 cc) , (4.7) 

where Tcs is the color symmetric piece 

Tcs{99 - cc) = d'^''^j^^e,{g,)e,{g2)i^ze^''"^{g^ - g2)aPp i'^T^V 
+ [2e'^'^"^((7i - ^72)ae*^'Aj + 4m,(Af A^^ + A^A^^) 

- — (^71- m2gi^k) - 2g^h^ + g^^{g^ - g,) ■ A,) ] g^^V^T^r^j + ■ ■ ■ , (4.8) 

and where Tqas is the color antisymmetric part 
TcAs{99 - cc) = ^r^'^^A_,^^g^)e^{^g^) 

2^;2^AJ - 2gm + g^^{g, - g^) ■ A^ ^^'T^V + " • • • (4-9) 
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The dots ■ ■ ■ are included as a reminder that we are writing only the first term in the 
infinite and non-convergent series in q^/ m^. This is justified since the factorization formalism 
requires only that we know the algebraic form of T for small q. 

We note that at leading order in the non-relativistic expansion, the amplitude for produc- 
ing cc{8,^Si)) in Figure 3b, Eq. ( [4.9[ ) cancels against the contribution from the gluon-fusion 
graphs in Figure 3a, leaving only the pieces in Eq. (|4.8|) . 

We now turn to the task of isolating the terms in Eq. ([4.8|) that contribute to the 
production of cc pairs in the states 8,^ So, 8, ^Pq and 8, ^Ps- Be ginning with 8, ""^^o, we note 
that terms associated with the angular-momentum quantum numbers ^5*0 will have Pauli- 
spinor bilinears of the form C,^t], q^^^T], etc. It is an easy matter to pick out the ^5*0 pieces 
in Eq. 



Tigg ^ cc(8, ^^o)) = ^d^^'^^jT^, e,{g,)e,{g2) e^'^^^ig, - g^^Pp i^T^^ (1 + ■ ■ • (4-10) 
Squaring Eq. ( |4.10| ), and performing the sum-average over spins and colors, we obtain 

\T{gg cc(8, 'SoW = ^^^'^V V^T^ (1 + ■ ■ ■) • (4.11) 
Inserting this into the expression for the cross section given in Eq. ( |3.12| ) we obtain 

a{gg ^ cc(8, ^^o)) ^ ^^^^ J m ' E{) ^T^v v'TX (1 + ■ ■ ■) , (4-12) 

which is the result in full QCD for the production rate of an on-shell cc(8, ^So)- 

The next step consists of calculating the same production rate in perturbative NRQCD. 
This is given by 

a{gg - cc(8, ^^o)) = / {0\Of'\'Som + 

J2 (0|x"^^>|c(s,"^,q)c(s,m, -q)) 

ssmfn 

,m,-q)|'*tT-x-|0> + .-- 

(2.)3 < f*T=',,'n+.... (4,13) 







(27r)3 


ml 


d^q 




(27r)3 




{c{s 


,m, q)c( 


d^q 
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Matching the result from full QCD in Eq. (|]T2]) to that from NRQCD in Eq. (|4T3| ), we 
obtain 

from which we conclude that 

a{gg ^ + X) = M^^^pfi^^ s,m 



5{s-Aml){Q\Ot{'S,m. (4.15) 



96m^ 

Next we isolate the individual P-wave contributions. This can be accomplished by first 
noting that any direct product of cartesian vectors may be written as 

a'V = ^6'^ a ■ b + ^e'^\a x h)'' + a^V^ , (4.16) 

where a^'V^ = {a'V + a^U) /2 - a • b 5^73. Thus we are able to decompose the factor q^^^a^r] 
appearing in the second term on the right- hand- side of Eq. (|4.8| ) into a scalar component 
that is identified with the ^Pq state, a vector component identified with the ^Pi state, and a 
symmetric-traceless tensor component identified with the ^P2 state. We wish to emphasize 
that the procedure of decomposing the amplitude into seperate angular-momentum configu- 
rations, ^5*0, ^5*1, ^Pq, ^Pi, etc. can be carried out to any order in the q^/wi^ expansion. It is 
always possible to decompose the amplitude for the production of a heavy quark-antiquark 
pair into pieces of the form 



r{^- 


- ccCSo)) 


r(.- 




T(z- 


- c-cCP^)) 


T(.- 




T(.- 


- c-cep,)) 


T(z- 





1* ^-t. 



m2 



P2)) = A'^ e^g'V^lr^ (l + a^ + --] 
\ mi J 
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where the factors A, A'', A^^ are independent of q. 
The ^Pq contribution can be written as 

r{gg ^ cc(8, 'Po)) = d^''-^ e,{g,)e,{g^) {-2g^''ml + P^g'() {^q ■ a + 



Squaring and performing the sum-average over spins and colors, we obtain 



(4.17) 



\n99^cc{8:'Po))\' = ^^^, eci-cT^V V^^.aTX (1 + ■ ■ ■) ■ (4-18) 



Inserting Eq. ( |4.18| ) into Eq. ( |3.12| ) we obtain 

a{gg - cc(8, ^Po)) = ^^^^ / ^ ' " ^T^V r/^q " ^T^^ (1 + " " ■) , 

(4.19) 



which is the result from full QCD. As to the perturbative NRQCD result, this is 



'j{gg^cc{8,^Po)) = J 



l^s^^ E ^(0|x^(-^B-a)T>|c(.,m,q)c(.-,m,-q)) 



(27r)3 


4 ^ 


d^q 


Fs^Po) 


(27r)3 




(c(s 


m, q)c(s 


d^q 


Fsi'Po) 



I 
2 



(27r 

Matching Eqs. (|4T9| ) and (^I20D , we obtain 



^8('Po) = ^5(3-4m^) (4.21) 



from which we conclude that 



-^8(^-Po) /nl/o^^3, 

6is-Aml){0\OtCPom. (4.22) 



a(<;^ - + X) = ^V^(0|03^(^Po)|0) 



(2me)5 

The 3Pi piece in Eq. ( ^.8| ) is found to vanish exactly. 
Extracting the ^P2 contribution from Eq. (^4.8|) , one obtains 
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X <^ 2m,Af a; - 



{91 ■ A)i 



rrir 



(gi^^';-g'i^^^+g^''{g,■^), 



x^^gMT^r^ +•••. (4.23) 
Anticipating that the cross section will involve an integral over (i^q, we note the identity 

I 7^ ^^Q^'^'^T^V V^q^'<y'^TX , (4.24) 



(27r) 



5 I 



2 3 ; J (27 ^ 

which, upon squaring Eq. (|4.23| ) and performing the sum-average over spins and colors allows 



us to write 



3 o m2 



(2vr) 



/ (Ct^ ^^^''^''^'^ r^^glV^lT-^e (1 + ■ ■ ■) • (4.25) 



24(2me)6 J (27 

Substituting Eq. (|]25|) into Eq. ( pl2D , we obtain the full-QCD result 

^(^7^7 - cc(8, ^P^)) = ^ / ^(^/ - E^) ^q^'cr^^T^V V^q^'a^^^^ (1 + " " ") • (4-26) 

On the other hand, the expression for a{gg — > cc{8,^P2)) calculated to lowest order in 
perturbative NRQCD is 



aigg^cc{8,'P2)) 







(27r)3 








(27r)3 
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{c{s 


,m, q)c( 


d^q 


Fs{'P2) 



{o\of''\'P2m + 



E (0|X^ (-^ B ' j T>|c(s, m, q)c(s, m, -q)) 



■2^ 



+ ■■ 



(4.27) 



Matching Eqs. ( |4.2(j| ) and ( ^.27] ) we obtain 



i^8('P2) 



24m^ 



S{s — Ami 



(4.28) 



from which we conclude that 
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= 3(2^' ~ ^""'^ (0|(^8n'i'2)|0) . (4.29) 
We can now write the proton- antiproton ■^/'-production cross section by convoluting the 
subprocess cross sections with the parton distribution functions. From the subprocess cross- 
section for qq —>■ ip + X , given in Eq. (|3.21|) , we have 



a{pp -^tp + X) = J dXifqlp{Xi) j dX2fq/p{X2) 



9^2 3 

^6{x,X2S-4ml){OtCS,)) + q^q, (4.30) 

where s = (Pp + Pp)'^ is the center-of-mass energy squared of the colhding proton-antiproton 
system. From the subprocess cross-sections for gg ^ + X , given in Eq. ( [4.15|) , Eq. ( [4. 221) , 
and Eq. ( [4.29| ), we have 



a{pp + X)= J dxifg/p{xi) j dx2fg/p{x2) 

(4.31) 

This expression can be further simphfied using the relation 

{Oti'Pj)) ^ (2J + l){Ot{'P,)) , (4.32) 

which holds to within corrections of relative order v"^. 

Finally we note a means of checking our results for the subprocess cross-sections gg — >■ 



cc(8, ^'^+^Lj) given in Eqs. ( [4.15|) , (f4.22|) and ( |4.29| ). This check entails first noting the 



surprising fact that the amplitudes for gg cc{8,'^^'^^Lj)) are proportional to those for 
gg — > cc{l,'^^~^^Lj)). Thus, our results for octet production rates should be related — by 
an overall color-factor and a replacement of matrix elements - to results computed from 
the same subprocesses for rjc and Xc production in the color-singlet-wavefunction model. 
These color-singlet results can be obtained from our color-octet results by first performing 
the replacements 
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{OteP2)) - f (Of . (4.33) 

To convert these expressions into the language of the color-singlet-wavefunction model, one 
makes the further replacements 

Zn 

(OrePj)) - ^^^I±^ |i?^(0)P . (4.34) 

One observes that, indeed, these operations yield results which agree with those presented 
in Ref B. 



V. CONCLUSION 

In this paper we have presented a calculation of the hadronic ip production cross-section, 
carried out within the framework of the NRQCD factorization formalism of Bodwin, Braaten, 
and Lepage. We have explicitly shown how to put into practice the "matching" procedure, 
which allows a determination of the short-distance coefficients appearing in the factorization 
formula. By remaining loyal to, and expanding upon, the matching program briefly described 
in Ref [|l| , we have revealed the true spirit of the BBL formalism for quarkonium production. 
We have also demonstrated how, operationally, one may systematically include relativistic 
corrections. 

We have obtained the following results for the leading order subprocess cross sections 

due to the color-octet mechanism: 

2^2 3 

^ilQ ^ ^ S{s ~ 4ml) {0\OtCS,m 

^(99 - ^ S{s - Ami) {{^\Ot{'S,m + ^(0|0|'(^Po)|0) + ^.{^P^'P^M^ ■ 

As a check on these results, we have transformed the above expressions into color-singlet- 
wavefunction results for r/c and Xc production, and have found them to be consistent with 
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previous work. We wish to point out that our expression for a{gg ip) differs, by a factor 
of 3, from an analogous expressions obtained in Ref [0]. 

The formula summarized above are of great relevance in the parametrization of quarko- 
nium production, and will be useful in future analyses. In particular, these expression, along 
with the previously calculated color-singlet result, appear in a leading order calculation of 
hadronic ip production at fixed target experiments. In fact there are discrepancies between 
theoretical predictions and experimental results for ip production in pion-nucleon interac- 
tions 1^ that will be resolved by including, in the theoretical predictions, contributions from 
color-octet subprocesses 0. 

Also note that the NRQCD matrix elements (0|C^(3^i)|0), (0|C»^(i5o)|0), (0|0^(3Po)|0), 
and (O|0^(3P2)|O) that appear in our calculation, in the linear combinations summarized 
in the above formulae, will appear in different linear combinations in a calculation of ip 



production in hadronic colliders |10 



It must be pointed out that our result for a{gg —>■ ■p) is readily converted into an 
expression for the forward photoproduction rate a{'yg 'p), which is the subject of a work 
in progress . 

Finally, we mention that while this work was in progress, there appeared an article by 
Peter Cho and Adam Leibovich in which very similar work was presented. These authors 
calculated the leading order -p production subprocess cross-sections that we present here, 
and we agree with their results. We must point out however that the approach to matching 



in [10 1 is quite different from ours, and that the two articles are complementary. 
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TABLES 

TABLE I. The t '-scaling of the NRQCD color-octet matrix elements, relative to S-wavo baseline. 

(0|Of ((25+i)Lj)|0) Scaling 

{o\otepjm/mi 

TABLE II. The scaling with as{mc) of the short-distance coefficients 



i + j ^ cc{n) + X Scaling 

99 cc(i, ^Si) + g asinicY 

qq cc(8, ^^i) "^(mc)^ 

99 cc(8, ^So,^Po,2) as{mc)^ 
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FIGURES 
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(b) 

FIG. 3. Color-octet subprocesses gg cc(8, ^5o, ^Si, '^Po,2)- Figure 3a: Gluon fusion 
Feynman diagram, which produces cc(8, ^^o, '^Si, '^Po,2)- Figure 3b: Feynman diagram for 
gg ^ g ^ cc(8, "^Si). At leading order in the non-relativistic expansion, the amplitude for cc(8, ^Si) 
production in Figure 3a cancels against that in Figure 3b. 
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